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Abstract: A fully classical explanation of the nonhydrogenic ionization threshold for low angular 
' momentum Rydberg states of Alkali-metal atoms in a linearly polarized low frequency monochro- 

\ matic microwave field is given: the classical equivalent to the quantum rate-limiting step, which is 

^S) ■ responsible for the scaling and which according to the literature initiates what then continues 

^ ' as essentially classical diffusion, is found. 
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I. INTRODUCTION 



Laboratory experiments [1,2] have shown that while the ionization threshold of high angular momentum Rydberg 
^ Alkali-metal atoms in a linearly polarized monochromatic microwave field with intensity F and frequency lj <C n^'^ , 
Q where uq is the initial principal quantum number, follows the hydrogenic behaviour Ftu — (IGj^q)"^' the threshold 
for low angular momentum states behaves as (Sug)"^. The quantum numerical simulations of A. Buchleitner and A. 
Krug [3-5] have substantially confirmed this picture, even though they were not able to definitively rule out for low 
angular momentum states a Fth ~ (OOtiq)"^ behaviour [5] for which on the other hand no theoretical explanation is 
known. 

In the literature it is often found the opinion that this (Sng)"^ threshold behaviour is nonclassical: see e.g. Ref. [6] 
', where the process is described as an essentially classical diffusion one initiated by a rate- limiting step which is a Landau- 
rS i' Zener transition at a quantum avoided crossing. On the other hand we are in a regime where classical behaviour can 
be expected on the relatively short time scales of the experiments: ionization of high angular momentum states follows 
the classical hydrogenic behaviour Fth — (16?T-o)~"'^i dynamical quantum localization does has no relevance at such low 
^ \ frequencies as the quantum delocalization border is below the experimental ionization threshold [7,8]; tunneling [9] 
\^ 1 and multiphoton processes happen at much longer time scales; and finally avoided crossings -the quantum resonances 
' which are the foundation of the quantum explanation given for the (Stiq)"^ threshold- are in most cases clearly related 
. to easily identifiable classical resonances, the magnitude of the splitting being proportional to the width of the classical 
resonance zone [10]. The aim of the present paper is to show that a fully classical explanation of the (Stiq)"^ threshold 
exists and thus add to the evidence in favour of such a scaling: on the energy surface of a Rydberg Alkali-metal atom 
in a static electric field the atomic core potential induces resonances; their overlap -which we shall show happens at 
a field intensity F k, (3ng)~^- brings chaos [11] and if the system is chaotic, a slowly varying field cannot be followed 
Q ' adiabatically [12,13] by the system and energy diffusion will ensue. 
• i-H ' The paper is thus organized: section II summarizes the quantum explanation found in the literature; in section III 
classical numerical simulations are presented which suggest that a classical explanation is also possible; in section IV 
the numerical results of section III are explained through a study of the phase space structures on the constant energy 
Q_i surfaces of an alkali-metal atom in a static electric field. Finally in section V the results of section IV are applied to 
the case of a slowly oscillating harmonic field to derive the experimentally observed ionization threshold. Section VI 
summarizes the results obtained. 



II. THE QUANTUM PICTURE 



The standard quantum explanation of the (3nQ)^^ behaviour is given [6,14] in terms of core-induced Landau-Zener 
interactions [15] between states of Stark manifolds with principal quantum number differing by 1 and field-induced 
Demkov-like interactions [16] within each manifold, as exemplified in the qualitative plot shown in Figure 1, which 
shows the optimal ionization path for a rig 3> 1 s-state as the electric field slowly oscillates at the threshold amplitude 
F = (3ng)~^. Interactions which can be approximated by the Landau-Zener model (characterized by field dependent 
diagonal matrix elements and constant off-diagonal ones) are marked by circles; interactions which instead can be 
approximated by a Demkov-like model (characterized by constant diagonal matrix elements and field dependent 
off-diagonal ones) are marked by squares. As the field rises the no,s state first undergoes an avoided crossing of 
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the Landau-Zcncr type with the highest energy state of the (no — 1),to = manifold and then, at _F = (Stt-q)"-'^, 
with the lowest energy state of the no,m — manifold; when the field decreases the population distributes among 
these three states and -as the field approaches zero- among all the states of the two manifolds, due to Dcmkov-like 
interactions. When the field rises again, the peak field is not high enough to couple the (no — l),m = manifold to 
the (no — 2), m = one; on the other hand several states of the no, m = get to be coupled to the (no + 1), m = 
manifold, so that diffusion toward higher energies and eventually ionization ensues. The actual ionization rate is not 
easily evaluated as the population fraction exchanged at each avoided crossing depends in a nontrivial manner on the 
oscillation frequency of the field and its peak intensity [17]. 

Since, in the experiments performed up to now, high angular momentum states arc not excited selectively in the 
quantum azimuthal number m, only a small fraction of their population is initially in the m = 0, ±1 Stark manifolds, 
which are those displaying significant core- induced avoided crossings [6], while most of the population is in the 
high m manifolds where the avoided crossings arc too narrow to induce nonhydrogenic behaviours. Manifolds with 
the same principal quantum number n but different m are not coupled by the electric field; the population of high 
angular momentum states is therefore trapped in the high m manifolds and the states themselves display a hydrogenic 
ionization threshold Fth — (16no)~^. 



III. PRELIMINARY CLASSICAL EXPLORATION 



As we have seen in the quantum description above, because the electric field varies very slowly, most of the evolution 
is adiabatic, and the results obtained for Alkali-metal atoms in a static electric field can be used to analyze the case 
of a slowly varying field. Our first step will therefore be a classical study of the static field case. Prom the discussion 
in the above section it is also clear that the case m = is exactly the case where the nonhydrogenic behaviour we 
want to investigate appears most clearly; I have therefore chosen to restrict my study to such a case. 

The motion of the perturbed system being restricted to a plane, I consider a two dimensional Alkali-metal atom 
model [18] in the {a;, z} plane: 

Ho = l{pl+pl)-l-^, (1) 



where px and pz are the conjugate momenta to the spatial coordinates x and z. The first two terms in eq. (1) represent 
the Hamiltonian for a Hydrogen atom with zero z component of the angular momentum = (that is: azimuthal 
quantum number m = 0) and the last term is the simplest known model for the nonhydrogenic core potential, /3 being 
the core charge and a a parameter proportional to the inverse of the core radius [19]. The values of the parameters 
which give the best fit to lithium have been found to be a = 2.13 and (3 = 2 [18]. 

We now add to the free atom Hamiltonian (1) a static field potential V = Fz; the system described by the 
Hamiltonian H = Hq + V is invariant under the scaling x ^ x/uq, z ^ z/hq, Px —^Px^ajPz Pzn-o, H Huq, 
a — > ariQ, 0^/3. Such a scaling would result in a hydrogenic threshold behaviour as can 
be easily seen from the scaling for the field intensity F: on the other hand the core parameter a is fixed for each 
atom type and therefore cannot be changed. Changing no without changing a means changing the relative dimension 
of the core to the orbit average radius and results in a non-hydrogenic threshold behaviour as can be seen from 
the following preliminary numerical results, presented -in accordance with the current literature- in semiparabolic 
coordinates scaled to the initial quantum number no: 

Vr + z 

u = (2) 

no 



^/r- z 

V = — , 3 

no 

and in their canonically conjugate momenta p„ and p-^. 

My findings are summarized in Figs. (2) and (3) where single orbit Poincare surfaces of section (SOS) for no = 40 
and no = 320 are shown. A SOS for a system orbit is given by the points where the orbit crosses a given plane with 
positive (negative) velocity. Those shown are for the {v,py) plane at u = with positive velocity u = du/dt > 0. 
A comparison of the right side of each figure, for which Fn^ = 0.04, with the left side where i^no = 0.32 clearly 
illustrates that the fraction of the energy surface a single orbit is able to explore varies as Fuq; in particular it can 
explore the entire surface when F = (3no)~^, as can be seen from the left hand side of both figures. The form taken 
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by this exploration is on the other hand different in the two cases: outside of the central region, where the chaotic 
behaviour first becomes evident [18], the no = 320 orbit appears much more regular than the no = 40 one; this is a 
consequence of the fact that far from the nucleus the system is hydrogen- like and therefore scales as Fnff. 

Our aim is to now understand why this happens; we shall do it by applying Chirikov's resonance overlap criterion 
for transition to chaos [11] in a form suitable to autonomous systems [20]. 



IV. THE CONSTANT FIELD CLASSICAL MODEL 



The quantum mechanical studies quoted in section II have shown that the relevant feature in the quasistatic 
ionization of alkali-metal atoms is the core-induced avoided crossing between the highest lying state of the no Stark 
manifold with the lowest one of the no + 1 one (in the hydrogen atom, the levels cross because of symmetry [21]). 
This avoided crossing takes place at an electric field strong enough that the electric field potential dominates over 
the core potential; it is therefore convenient to consider an (hydrogen-like) atom without core potential in a static 
electric field as our unperturbed system and the core itself as our perturbation; this puts us out of the range of the 
weak electric field approximation used in Refs. [22,23]. 

To apply Chirikov's criterion it is convenient to write the Hamiltonian in action-angle variables, so that we can 
explicitly calculate the characteristic frequencies of the unperturbed system. To do this we first write the Hamiltonian 
of the unperturbed system (a two-dimensional hydrogen atom in a static electric field) 

Hh = 1{pI+pI)- , ^ +Fz (4) 



in parabolic coordinates defined by 



(5) 

(0) 

and their conjugate momenta and p,,. The resulting Hamiltonian 

separates [12] and is therefore regular (non-chaotic). We can now pass to the action-angle variables {/, /i, A, /i} defined 
by the equations 

e = 2//i(l-sinxi) (8) 
77 = 2/(/-/i)(l-sinx2) (9) 
1 cosxi 

^« = 27T3iI^^ (10) 
,, = J-.J^, (11) 

27 1 - sm X2 

where 1\ < I, and xi and X2 are auxiliary angles defined by 

h I-h , TT 

A = -yCosxi ^cosx2-X2+2 (12) 

M = X2-Xi- (13) 

The two actions I and I\ are the classical equivalent of the parabolic quantum numbers n and n\ used in quantum 
mechanics to describe the interaction of a hydrogen atom with a static electric field [24] and A and /U are their 
respective canonical angles. 

In these coordinates, the Hamiltonian (7) reads, in first (linear) approximation, [12] 

H = -^^+\FI{2h-I). (14) 
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Wc have considered only the hnear term of the electric potential as the relevant crossing in the corresponding 
quantum system takes place in the linear regime of the Stark potential not only for the high angular momentum 
states but also for the low angular momentum ones which, due to the quantum defect, display at low electric fields 
only a quadratic Stark shift [25,26]. 

We can now add the core potential to the above Hamiltonian; to introduce it as a perturbation, we take its Fourier 
expansion: 

V = Ao,o(/, h) + T,k,k,>o2AkM (I, h) cos (fcA + fci/x). (15) 

where the coefficients Ak^ki{I,Ii) are the seniiclassical matrix elements of the perturbation [24]. For the resonant 
terms wc are interested in, we have -as we shall shortly see- fc <C /ci and k <IC aP; as long as /i is not too close to 
or / (the condition Ji, (J — /i) > k\/{2al) must be verified), we can therefore write these coefficients as: 



(16) 



■^xp[- — 1-7- rT] = Ta ' (17) 



27r/4a/i(l - /i) 2q!/i(1-/i)' /4 

where the IjiuYs are modified Bessel functions and we have introduced the scaled action Ji = Ii/I and the scaled 
index k\ = —k\/I. The approximations made to derive eq. (17) are discussed in the appendix. 

The scaled coefficients A^, _fej(/i) appear therefore almost independent from / and k. Figure (4) shows the depen- 
dence on kx of the exact expression for A^ eq. (16) for several values of I\ for lithium (a = 2.13 and /3 = 2): 
the matrix elements are significant only for ki ~< 1. 

Figure (5) shows instead the dependence of Ai._i(/i) on /i, again using the xact expression eq. (16): there is a 
wide peak around the central value I\ =1/2 where from eq. (17) we have 

-2 

_i(l/2) ~ —e— =~ 0.116, (18) 
' 7ra 

and drops to zero for Ji = and /i = 1. 

Resonances induced by the core will take place when the phase of one of the terms of the Fourier expansion (15) is 
stationary; this happens when the two unperturbed frequencies of motio for A and fjb, 

dH I , , 

io<,^ — = j^+^F{h-I) (19) 

0^1 = ^ = 3FI, (20) 

ail 



satisfy the resonance condition 



kiOQ + kioji = 0, (21) 



with k and ki two integers. 

Since, for F below the hydrogen ionization threshold, we have wq ^ <^i, the most important resonances (those 
having large stationary terms) will be for fc = ±1, so that fci is not too big; in particular from Figure (4) we see that 
the matrix elements Ai^-ki{I,Ii) are significantly large only for ~< /. Substituting k = 1 and imposing —ki < I 
in the resonance condition eq. (21), we obtain 

^ ^ 3P [i + lii-ii)y 

For the high values of the principal action / we are considering the second term in the denominator is negligible; eq. 
(22) therefore becomes independent of Ji and, since / is the classical analogue of the principal quantum number no, 
corresponds to the experimental ionization threshold F — (Stiq)^^. 

Obviously, the above argument only tells us that classical resonances are noticeably big only for F >^ (Sng)""^; 
whether they do overlap and thus generate global chaos on the energy surface is still an open question. To answer 
this question we have first to find the positions of the resonances on the energy surface. Let us therefore place ourself 
on a constant energy curve in {7, /i}: 
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Eo = ~^ + '\FI{2h-I)- 



(23) 



introducing the average action Jq = 1/\/—2Eq (corresponding to h = 1/2), and the scaled quantities / = I / Iq and 
Fq = FIq, eq. (23) now reads 



3Fo7''(27i - 1) + P - 1 = 0. 

In the same variables the resonance condition cq. (21) for k = 1 roads: 

3FoP(-fci + 1 - Ji) = 1. 
If we now eliminate I between eq. (24) and eq. (25), we obtain the condition 

3Fo(-fci + 2 - 3hf = (-fci + 1 - Ji). 

which, solved for ki, gives us 



-kf 



1 - 6Fo(2 - 3Ji) ± - 12fb(l - 2Ji) 
QFo 



(24) 
(25) 
(26) 

(27) 



Only the positive root is of interest to us and for 12i^o <C 1 we can expand the square root obtaining 
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3 - 3Fo(l - 2/i)^ + ~ +5/1-3 (28) 

61<o 

Since 7i e [0, 1], we have five resonances (six if l/(3Fo) is an integer) at equally spaced values of Ji between and 



1: 



h = 



3Fo 



modi 



J=(0),l,...5. 



(29) 



where the value of the index J has been put in parenthesis to indicate it is possible only for l/(3i^o) an integer. 
The corresponding indices fci will read 



1 



3Fn 



J-3 



where [.] denotes integer part. The resonance values for / will instead be 

1 



{l + 3Fo[|(j-,iLJ-2]} 



1/4' 



(30) 



(31) 



which, since we are considering SFq ~ I/Iq <C 1, we can approximate with 1. 

From Figure (5) it is clear that the widths of the first and last resonances will be zero (for l/(37b) an integer) or 
close to zero, while the other four will be of comparable widths. We now want to compare the average width of these 
resonances to the separation of two consecutive resonances, if the former is larger than the latter, the resonances 
overlap and we have chaos on the energy surface (Chirikov's resonance overlap criterion [11]). 

To apply Chirikov's resonance overlap criterion for transition to chaos to autonomous systems, we have to take 
into account energy conservation [20]: the resonance width to be considered is the one along the constant energy 
curve in containing the resonance center itself; from eq. (14) the versor f tangent to the energy curve reads 

f = — wo/l'^l} where cjq and ui are given by eqs. (19,20) and \u)\ = y/ul+ujf. Following Ref. [20], the 

resonant Hamiltonian in the restricted phase space of the tangent action J = [(/ — I^^^)uJi — {h — Ii'^)(Jo]/\i^\ and its 
conjugate angle tp = (Awi — fjuo)/\io\ where /^''^ and I^"^ are the actions at the {k,ki) resonance, reads: 



Hr = -aJ^ + 2Ak,k^ (/^'"\ 4"') cos (il^), 

fl = k\Lu\/u}i, 
a = far, 



(32) 

(33) 
(34) 
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where a is the matrix of the second derivatives of H with respect to / and /i. As |a;| ~ ujq, we have f2 ~ |fci| and |A;i| 
is the number of islands in the chain along the variable 
The half width of the (fc, fci) resonance in J is therefore: 



and since |a;| — wq, it is also the width in /i, while the width in / is negligible. 
Since i^o ^ 1 we have, dropping the index (r) to avoid too cumbersome a notation, 

/3F/\7-5 h\ fSFoP\' -5 + 3foP(l-2/i) .{WV 

' = [-wj [i^^'^- j = [-ir) [i+3Fom-iF+(3w " y-ir) ■ ^''^ 

Taking for ki the approximate average value from eq. (28) fci = — l/(3i^o) and A; = 1, so that fc — fci ~ — fci; taking 
also Ji = 1/2 (this latter being the value for which Ak^ki {I, h) is maximum, it will give an overestimate of the width) 
we obtain, using eq. (17) and remembering that for Ji = //2 we have / = 1, 



^,-^(/o,/o/2)^^e^. (37) 



The maximum resonance half width on the energy surface, scaled to the average principal action Iq, will therefore be 



W 2/2/3 ^ 

lo 3Fo/o V 57ra ^ ' 

Figure (6) shows W/Iq as a function of 3Fo/o for the lithium parameters. The dotted curve is eq. (38); the full curves 
are instead the exact scaled widths calculated at Iq = 10 for (from top to bottom) Ji = 0.50, 0.20, 0.10, and 0.05. The 
first overlap will happen when the Ii = 0.5 width becomes larger than 0.10, that is for 3Fo/o — 0.4; to have global 
chaos on the energy surface we instead need the extreme resonances to both overlap with their nearest neighbour and 
to reach the extremes in Ii of the energy manifold; this happens when the /i = 0.20 width becomes larger than 0.20, 
that is at 3Fo/o — 0.7, in reasonable agreement with the numerical results from section III: the deviation is within 
the usual factor 2 expected for Chirikov's criterion. Worth noting is the fact that, as the perturbation is the atomic 
core, the resonance width depends on the electric field F only indirectly, through the resonance index ki (see eq. 
(30)) and the energy curve shape. The first dependence gives the growth of the resonance width, the second one its 
slow decrease at high F where on the other hand the linear approximation for the electric field potential used in the 
Hamiltonian (14) is no more valid. 

While eq. (38) shows that the average behaviour depends on the product 3Fo/o only, the details depend on 3Fo 
and Iq separately: using eqs. (17), (29), (30), (35), and (36), we obtain, remembering that from eq. (31) I ~ lo, 



W 



lo 3Fo/o 



.3^1+-'- 

2T7\ — 



modi 



3Fo 



^'-^Lo.l'^'-'^^\mo.l\ (39) 



modi 



Fig. (7) shows the positions (eq. (29)) and widths (eq. (39)) of the resonances in Ji = Ii/I ~ Ii/Io as a function 
of 3Fo/o for a) Iq = 10 and b) Iq = 60: in both cases the general behaviour is the same: the first overlaps happen 
for Ji = 0.5 at 3Fo/o — 0.41 (see the detail in Fig. (8)) and the entire Ji range is covered by the resonances for 
3Fo/o — 0.7, but the vertical symmetry is much better for Iq = 60; this is easily explained by noting that, because 
of eq.(28), the speed of motion of a resonance toward higher 7i values with increasing 3Fo7o increases with Iq, and 
that in our approximation the resonance width depends from Ii only through ^^^^,^(7, 7i) and is therefore symmetric 
in 7i and 1 — 7i along the resonance curves {ki = const). 



V. THE VARIABLE FIELD CLASSICAL MODEL 



Let's consider a cyclical perturbation of a given system; an action is an (approximate) adiabatic invariant if the 
unperturbed frequency of motion of the associated angle variable is much larger than the fractional rate of change 
of the perturbation itself [12]; even though there is no change of energy over a whole cycle of the perturbation, the 
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energy of the system changes during the time the perturbation is varied. The opposite Hmit is that of sudden or 
diabatic variation: in this case the energy of the system does not change at aU during the time the perturbation is 
varied. 

In a few words: diabatic variation of a perturbation means that the motion of the unperturbed system is very slow 
compared to the rate of change of the perturbation and therefore a whole cycle of the perturbation brings no change 
to the system, which had no chance to evolve. Adiabatic variation means instead that the motion of the unperturbed 
system is very fast compared to the rate of change of the perturbation so that the system can follow the slowly varying 
orbit it is on and again a whole cycle of the perturbation brings no change to the system, apart from a vector phase, 
sum of a dynamical part (the integral over the pulse time of the istantaneous frequency for each action) and of a 
geometrical part (the so called Hannay angle [27]). 

Different orbits having different periods, adiabaticity depends on initial conditions. Two classes of orbits are 
intrinsecally non-adiabatic. The first one is that of the seperatrices: orbits with infinite period which separate regions 
of phase space characterized by different kinds of motion, e.g. rotation and libration for a pendulum. The other class 
is that of chaotic orbits, beacause of exponential instability. 

The really interesting case is the case in between the two above: in this case motion of the unperturbed system and 
rate of change of the perturbation are comparable and the system neither stays still nor evolves along the varying 
orbit it was initially on. It instead moves to different static perturbation orbits so that in general energy exchange 
will ensue. 

A particularly interesting such case is when the variation of the perturbation causes the system to non-adiabatically 
cross over (cither through a separatrix or a chaotic region) from a region of phase space to another where motion is of 
a different nature; in this case the phase acquired by the system during the adiabatic part of the evolution can play 
an essential function in determining the final energy (see e.g. [28]). 

In our case we slowly vary the electric field according to the sine law F{t) = F sin (oj/t) so that the fractional rate of 
change reads w// tan(a;/t); apart from the cases uift = Nn, N an integer, where it goes to zero, or ujft = {2N -\- l)7r/2 
where it goes to infinity [29], we can use the approximate value w/. 

Considering as our unperturbed system a Hydrogen atom in a static electric field the frequencies to be considered 
would be those associated with the two independent actions I\ and l2' in the electric field linear regime, they give 
the conditions 

< +3F7i = wo + wi (40) 

ujf^^-2,Fh=wo; (41) 

which, as in our case FI^Ii < FI^ <C 1, reduce in both cases to 

ujfl^ < 1. (42) 

The above condition is the exact definition of the regime we are interested in and is therefore always verified. 

In the case the unperturbed system is instead an Alkali-metal atom in a static electric field, we shall distingush 
two regimes: low electric fields F <C (Stt-q)"^; and threshold electric fields F ~ (3ng)^^, where, as chaos is global, 
adiabatic evolution is impossible everywhere. 

For low electric fields the situation is more complicated, as the phase space is mixed: deformed tori and resonance 
islands coexist in it; far from the resonance islands the two actions 7i and I2 are still good approximations of the 
actual actions of the system and we still have adiabatic evolution for uj/I^ ^ 1; in the regions on the border of the 
resonance islands the motion is instead chaotic and no adiabatic motion is possible [30]; finally, inside the islands 
the system is approximated by a pendulum one [11], from eq. (32) we therefore have that, close to the center of an 
island, the frequency of motion around the center of the island itself is Ui = fl^/2aAk^k^ i^jh), so that the adiabaticity 
condition reads: 

a;,7^«^Fo^; (43) 
using eq. (38) as an approximate evaluation of W this condition becomes: 



, 1 /10/3 -1 V3 -1 , 
/ V 7ra / 

where in the last expression we have introduced the parameters for Lithium. As 7 ^ 1, the condition is a very 
stringent one and is never verified for the range of parameters we are interested in; adiabaticity is therefore not 
possible also inside the resonance islands themselves. 
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Let us now see what happens when the electric field in our system is slowly increased: from Figs. (7) and (8) we 
see that the resonances appear at /i = and m,ove to higher values of Ii [31] growing in size up to Ii = 1/2 and then 
decrcasising again. This might at first sight suggest as a ionization mechanism the trapping of orbits in the resonance 
islands which would then transport them and finally release them either when the resonance islands themselves break 
because of overlap with nearby resonances [32] or when the island size decreases (either at high 7i for increasing field, 
or at low Ii for decreasing field). Such a mechanism could in principle work for F < (3no)~^; it would on the other 
hand require adiabatic (or quasi-adiabatic) evolution within the resonance islands, which as we have seen is far from 
the case for the range of parameters we are interested in: as the motion inside the resonance island is at its fastest 
still much slowlcr than the rate of field change, the islands themselves will not be seen by the orbits. 

The only possible mechanism is therefore diffusion in the chaotic sea, which happens only for F ~ (Stiq)"^. As it 
was the case in the quantum description summarized in Sect. II, an actual evaluation of the diffusion and ionization 
rates does not appear easy: the only data we have for comparison are the numerical ones in Ref. [5] and they show that 
over long interaction times the decay of the survival probability is neither exponential, as expected for a completely 
chaotic phase space [33], nor algebraic, as expected in mixed phase space [34]. On the other hand, for short interaction 
times, while for the near-threshold plot in Fig. 6.8 from Ref. [5] the figure resolution is too low to extract any data 
for comparison, the high field plots seem to indicate an initial behaviour still in accordance with the classical decay 
time to = ncoy^ / IqF'^ given in Ref. [35] for w/ > Iq^. 



VI. CONCLUSIONS 



We have shown here that while the ionization of excited hydrogen atoms in a quasistatic monochromatic field is 
regular, that of low angular momentum alkali-metal atoms is chaotic. Core induced resonances overlap at F ~ (3nQ)~^, 
thus triggering chaotic diffusion which results in ionization. The motion within the single resonance islands being 
extremely slow, the islands themselves are not seen by the system during the pulse and therefore allow no efficient 
path to ionization. 

The present paper concludes my study in terms of nonlinear dynamics of the ionization mechanisms for excited 
Alkali-metal atoms in microwave fields initiated in Ref. [8] and prompted by the individuation of three different 
regimes in the numerical studies by A. Buchlcitncr and A. Krug [3 5]: in the high frequency regime OJ/I^ > 1 the 
ionization threshold is determined by the same quantum localization mechanism as for excited Hydrogen atoms [35] ; 
in the intermediate frequency regime w//^ '~< 1 quantum localization again detemines the ionization threshold, but 
the Alkali- metal atoms quantum defects raise the dimensionality of the problem [8]; finally, as we have shown in the 
present paper, in the low frequency regime lj/I'^ <C 1, classical dynamics, which for excited Hydrogen atoms essentially 
determines the ionization threshold alredy for l^,'// ' < 1, determines the ionization threshold for excited Alkali- metal 
atoms too, which, due to core effects, is chaotic and much lower than the Hydrogen one. 
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APPENDIX A: EQUATIONS OF MOTION IN SEMIPARABOLIC COORDINATES 

We derive here the exact expression of the coefficients Ah,ki{I,h) of the Fourier expansion of the core potential 
and discuss the approximation eq. (17) given in the text: the coefficients Fourier expansion 



V ^ -^^-==- = ^kM^kM U> /i)e'(^^+'=^'') (Al) 
are given by the integral 

1 f'^^ r'^^ aVx^+z^ 

A,,,,{I,h) = -—^ dX ^ e-(''^+^^^). (A2) 

{■^T^r Jo Jo VX^ + 



8 



To evaluate it, we change integration variables from A and /U to Xi and X2'- from (12) and (13) the Jacobian determinant 
of the transformation reads 



D{xi,X2) = y smxi H ^smx2 - 1; 



(A3) 



using then eqs. (5), (6), (8), and (9), we obtain [36]: 

^k,ki {I, II ) = 

1 /3e-"^' '•2- 



(27r)2 P Jo Jo 



COSX2 + (fe-fel)X2 _ 



P \aP-k 



k/2 



4, (j^/aH^ - h-k, (^-^^/a^P - , (A4) 



where the Ij{yy& are modified Bessel functions. 

Because for integer indices the modified Bessel functions have the property I-j{y) = Ij{y), it is clear from eq. (A4) 
that A-k,-ki{I,Ii) = Ak^kiiI,Ii)- Eq- (Al) can therefore be written as: 



V = Ao,o(/, h) + Sfe,fe,>o2Afe,fei (/, h) cos {kX + kijj). 



(A5) 



which is the expression given in the text, eq. (15). 

If we now introduce the scaled action 7i = /i/7 in eq (A4), we can easily see that, if A; ^ aP, the arguments of 
the modified Bessel functions scale approximately as P . Moreover, whenever the index j is small with respect to the 
argument y, we can use the expansion [37] 



rhry 



if-1 {4f-l){4f-9) {4f-l)i4f-9){4f-25) ' 



8y 



2!(8y)2 



3!(8y)3 



which, if j is itself large, can be further simplified as 



rMy 



1 



2y 2!(2y)2 3!(2y)3, 



rMy 



(A6) 



(A7) 



Introducing the scaled index ki = —k\/I, assuming <C fci and 7i, (7 — 7i) > k\/{2al), and substituting eq. (A7) in 
eq. (A4), the matrix elements read 



(A8) 



which is the expression given in the text. From eq. (A8) we see that Aj, _f,_^{I,Ii) ^ 7""*; to give an idea of the 
residual dependence on 7, Figure (9) compares the scaled matrix element calculated using eq. (A8) with the exact 
ones for 7 = 10 and 7 = 20: even on a logarithmic scale the differences are small . 
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FIG. 1. (Color online) The optimal ionization path for a no S> 1 s-state in a quasistatic microwave field oscillating between 
iiFmax ~ (3no)~^. The rectangular box indicates the Demkov-like interaction region for the no manifold; circles mark instead 
the Landau-Zener interaction regions (to avoid overcrowding the picture, only the extreme relevant avoided crossings between 
the no and (no + 1) manifolds have been shown). 
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FIG. 2. The Poincaro SOS in v and p„ for no = 40 and Fn§ = 0.04 (right) and Fn§ = 0.32 (left). As the SOS is symmetric 
for V —V, only half of each SOS is shown 
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FIG. 5. The dependence of ^i,_i(7i) on h for a = 2.13 and /3 = 2. 
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FIG. 6. The dependence of W/Io on SFoIo for various values of Ii: from top to bottom Ji = 0.50, 0.20, 0.10, and 0.05. the 
dotted hne is the approximated width calculated using eq. (38). 
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FIG. 7. Position (eq. (29)) and width (eq. (39)) of the resonances in 7i = c:± Ii/Io as a function of SFqIq for a) lo = 10; 
b) 7o = 60. 
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FIG. 8. Detail of Fig. (7)b for the values of 3Fo7o where the resonances first overlap. 
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FIG. 9. (color online) Comparison of the exact dependence of Ai _j._^{Ii) on ki with the approximation using eq.(A8). a) 
7i = 0.05, full curve eq.(A7), dash: exact result for / = 10. b) Ji — 0.5, full curve: eq.(A7), dash: exact result for / = 10, dot: 
exact result for / — 20. the parameters are those for lithium: a = 2.13 and (3 = 2. The main figures are on a logaritmic scale, 
the insets on a linear one. 
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